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Abstract 

In this paper we first calculate the form factors of B ^ (vr, K) and Bg ^ K transitions 
by employing the perturbative QCD (pQCD) factorization approach with the inclusion of the 
next-to-leading-order(NLO) corrections, and then we calculate the branching ratios of the cor- 
responding semileptonic decays B/Bg — )• {n, K){l~^l~ ,liy,i>r') (here / denotes e,fj, and r). Based 
on the numerical calculations and phenomenological analysis, we found the following results: 
I (a) For B — t- {tt,K) and Bg ^ K transition form factors F^^j^^'i^q^), the NLO pQCD predic- 

I _ tions for the values of -Fo^+^t(0) and their g^-dependence agree well with those obtained from 

other methods; (b) For B^ — )• 'K^l~Di,K^l^l~ and B~ — )• 7T^l~i>i,K~l~^l~ decay modes, the 
• NLO pQCD predictions for their branching ratios agree very well with the measured values; (c) 

By comparing the pQCD predictions for Br{B^ — )• 7r+/~P;) with the measured decay rate we 
extract out the magnitude of Vub- \Vub\ = (3.80^Q5Q(t/ieor.)) x 10^^; (d) We also defined several 
ratios of the branching ratios, Ru,Rc and Rni,N2,N3, and presented the corresponding pQCD 
predictions, which will be tested by LHCb and the forthcoming Super-B experiments. 
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I. INTRODUCTION 



The semileptonic decays B — )■ {n, K){l^l^ Ju,!/!?) and Bg — )■ K{l^l^ ^W^vv) with 
/ = {e,fi,tau) are very interesting B/Bg decays modes and playing an important role in 
testing the standard model (SM) and in searching for the new physics (NP) beyond the 
SM, such as the determination of \Vub\ and the extractions of the transition form factors 
Po,+,'v{'f) of B/Bg meson to pion and/or kaon. For the charged current B/Bg — )■ PW 
decays, the "Tree" diagrams provide the leading order contribution. For the neutral 
current B/Bg — > Pl^l^ and Pui? decays, however, the leading order SM contributions 
come from the photon penguin, the Z penguin and the W^W~ box diagrams, as shown 
in Fig. 1, where the symbol © denotes the corresponding one-loop SM contributions. 

On the experiment side, some decay modes among the all considered B — )■ P{l^l~ , lu, vv) 
decays have been measured by CLEO, BaBar and Belle experiments The Bg — > 

K{l^l~ ^Iv^vv) decays are now under studying and will be measured by the LHCb and 
the forthcoming Super-B experiments 0, [^. 

On the theory side, the considered semileptonic decays strongly depend on the values 
and the shape of the B/Bg — )■ P form factors. At present, there are various approaches to 
calculate the B/Bg ^ {'^,K) transition form factors, such as the lattice QCD technique 
[8], the light cone QCD sum rules (LCSRs) liQ, as well as the perturbative QCD 
(pQCD) factorization approach |13l-ll7l|. The direct perturbative calculations of the one- 
gluon exchange diagram for the meson transition form factors suffer from the end- 
point singularities. Because of these end-point singularities, it was claimed that the 
B — 7- P transition form factors is not calculable perturbatively in QCD Jl8| . 

In the pQCD factorization approach jl9|, however, a form factor is generally written 
as the convolution of a hard amplitude with initial-state and final-state hadron distri- 
bution amplitudes. In fact, in the endpoint region the parton transverse momenta /ct 
is not negligible. If the large double logarithmic term ag\n^{kT) and large logarithms 
asln^(x) are resummed to all orders, the relevant Sudakov form factors from both the 
/ct resummation and the threshold resummation 20|, |2l| can cure the endpoint singularity 
which makes the perturbative calculation of the hard amplitudes infrared safe, and then 
the main contribution comes from the perturbative regions. 



In Refs. [13|, |16|, for example, the authors calculated the P — )■ tt, p [13[ and B S form 
factors 161] at the leading order by employing the pQCD factorization approach and found 
that the values of the corresponding form factors coming from the pQCD factorization 



approach agree well with those obtained by using other methods. In a recent paper [22 



Li, Shen and Wang calculated the next-to-leading order (NLO) corrections to the P — )■ vr 
transition form factors at leading twist in the factorization theorem. They found that 
the NLO corrections amount only up to 30% of the form factors at the large recoil region 
of the pion. 

In this paper, based on the assumption of the SU (3) flavor symmetry, we first extend 



the NLO results about the P — )■ vr form factors as presented in Ref. [22[ to the cases of 
B ^ K and Bg ^ K directly, and then calculate the g^-dependence of the differential 
decay rates and the branching ratios of the considered B/Bg semileptonic decay modes, 
and furthermore extract \Vub\ based on our calculations. 

This paper is organized as follows. In Sec. II, we collect the distribution amplitudes 
of the B/Bg mesons and the tt,K mesons being used in the calculation and give the 
fcT-dependent NLO expressions of the corresponding transition form factors. In Sec. Ill, 
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FIG. 1. The typical Feynman diagrams for the semileptonic decays B/Bg — t- 

{tt, K){l^l~ jW, i/jy), where the symbol in (b) and (e) denotes the flavor-changing neutral 
current vertex with V = j and/or Z boson. 



based on the /cx factorization formulism, we calculate and present the expressions for the 
B/Bg TT,K transition form factors in the large recoil regions. The numerical results 
and relevant discussions are given in Sec. IV. And Sec. V contains the conclusions and a 
short summary. 



II. THE THEORETICAL FRAMEWORK AND NLO CORRECTIONS 

For the sake of simplicity, we use B denotes both the B and Bg meson and P denotes 
final meson n or K from now on. As usual, we treat B meson as a heavy-light system. 
In the B meson rest frame, with the ttlb stands for the mass of B meson, we define the 
B meson momentum pi and the final meson P (say tt ot K meson) momentum p2 in the 
light-cone coordinates: 

Pi = -|(1,1,0t), P2 = ^v{0A,0t), (1) 

with the energy fraction rj = 1 — /m^ carried by the final meson (here q = Pi — P2) ■ The 
light spectator momenta ki in the B meson and ^2 in the final meson are parameterized 
as 

h = {xi^,0,kiT), k2 = {0,X2V^,k2T)- (2) 

Because of the final pseudoscalar meson moving along the minus direction with large 
momentum, the plus component of its partons momentum should be very small, so it's 
dropped in the expression of k2- But the four components of ki should be of the same 
order, i.e. 0(A), with A = ttlb — rrib, rub being the b quark mass. However, since ^2 is 
mainly in the minus direction with A;^ ~ 0{mB)-, the hard amplitudes will not depend on 
the minus component ki as explained below. This is the reason why we do not give k^ in 
Eq.([2]) explicitly. 

In Ref. the authors derived the /cx-dependent NLO hard kernel for the B ^ n 
transition form factor. We here use their results directly for S — )• vr transition processes. 
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and extend the expressions of -B — t- tt form factors to the ones for both B ^ K and 
Bg ^ K transitions, under the assumption of SU{3) flavor symmetry. As given in Eq.(56) 
of Ref. 22|, the NLO hard kernel H^^'^ can be written as 



i/(i) = F(xi,a;2,77,/i/,/i,Ci)i^^°^ 



a 



An 

,2 



21 



m 



13 



— In^- In-f + — ln^ + — ln'(xiX2) 



V Ci 

1 1 
+ - In xi H — In xi In xo + ( 2 In 



7 



m 



B 



16 



m 



B 



7 



Inxi 



+ 



- In r? 

8 ' 2 



lnx2 



15 7 

m r? m r? 

4 16 " ' 



m 



^In^^ ( 31n:^ + 2 1 + ^tt^ + ^ 



m 



B 



101 



219 



48 



16 



(3) 



where Ci = 25mB|22|, /i/ is the factorization scale and set to the hard scales ti or t2 as 
defined in the Appendix, r] = 1 — [pi — P2Y / is the energy fraction carried by the final 



meson, and finally the renormalization scale /i is defined as [22 



with the coefficients 
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In this paper, we use the same distribution amplitudes for B/Bg meson and for the vr 
and K meson as those used in Refs. [l^, 22-24 . 



')b{x, h) = Nbx'^{1 — xY exp 



and 



XTUb 



xrriB^ 

^Bs 



(6) 



(7) 



where the normalization factors A^'^^^j are related to the decay constants fs^^^ through 



dx4>B^^^ {x, 6 = 0) 



fB(s) 

2^6 ■ 



Here the shape parameter Ub has been fixed at 0.40 GeV by using the rich experimental 
data on the B mesons with Jb = 0.21 GeV. Correspondingly, the normalization constant 
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Nb is 101.4. For Bs meson, we adopt the shape parameter = 0.50 GeV with fs^ = 
0.23 GeV, then the corresponding normahzation constant is Nb^ = 63.67. In order to 
analyze the uncertainties of theoretical predictions induced by the inputs, we can vary the 
shape parameters Ub and ub^ by 10%, i.e., Ub = 0.40±0.04 GeV and lob^ = 0.50±0.05 GeV, 
respectively. 

For the vr and K mesons, we adopt the same set of distribution amplitudes 4>f{x) (the 
leading twist-2 ) and 0f'"^(x) with i = (tCjK) as defined in Refs. [13, [H, 25|): 



fix) = ^ x(l - x) \l + aiCf/'(t) + a^Cl^^t) + a^Cl^\t) 
v6 L 



2^6 

fi 
2^6 



1 + ( 30r^3 - IpI) Cl^'it) - 3 Ucos + ^p^(l + Qa^) [ C'/'it) 



x{l — x) 



X 7 3 



(9) 
,(10) 
(11) 



where t = 2x — l, p,r(x) = T^TT{K)/'m'o^^'' are the mass ratios ( here = m^/(m„ H-m^) = 
1.4 ± 0.1 GeV and = m\/{'ms + nid) = 1.6 ±0.1 GeV are the chiral mass of pion and 
kaon), a^' are the Gegenbauer moments, while C^{t) are the Gegenbauer polynomials 



Cfit) =3t, C',/\t) = \{3t^-l), Cnt) = ^ (5t^ - 1) , 
C'J^t) = i(3-30t2 + 35t^), Cf{t) = I^(l-I4t2 + 2U4) . 



(12) 



The Gegenbauer moments appeared in Eqs. (iQlfTTl) are the following [10|, |24 



0, af = 0.06 ±0.03, Og'^ = 0.25 ± 0.15, 



-0.015, 



TT,K 
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0.015, 
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(13) 



III. FORM FACTORS AND SEMILEPTONIC DECAYS 



A. Tr,K form factors 



The form factors for B ^ P transition are defined by [2 



m\ — m 



m 



B 



nil 



(14) 



where q = pi — p2 is the momentum transfer to the lepton pairs. In order to cancel the 
poles at = 0, -F+(0) should be equal to -fo(O). For the sake of the calculation, it is 
convenient to define the auxiliary form factors fi{q^) and f2{q^) 



{P{p2)\mi,qmB{pi)) = h{q')pi, + f2{q')p2. 



(15) 
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in terms of fi{q^) and f2{(f)-, the form factors -F+(g^) and -Fo('?^) are defined as 



1 + 



2 2 

— rrip 



+ lf2{q') 



2 2 

— mp 



(16) 



As for the tensor operator, there's only one independent form factor, which is important 
for the semi-leptonic decay 



{P{p2)\b{0)af,^q{0)\B{pi)) = i[p2^q^ - q^,P2u 



2FT{q') 
rriB + mp ' 



(P(p2)IK0)^,..75g(0)|fi(pi)) = e^,,^p?g 



TUB + mp 



;i7) 



The above form factors are dominated by a single gluon exchange in the lowest order 
and in the lar ge r ecoil regions. The factorization formula for the B ^ P form factors is 



written as [13|, [15 



{P{P2)\ m^,qmB{p,))=glCpKj dx,dx2Sk,^d'k2^^^'^^^Lf^ 

xe-^^-^(P(p,)|g;(l/)g;3(0)|0)e^'=-^(0|6.(0)gK^)|i?(pi))T/;"'- (18) 

In the hard-scattering kernel, the transverse momentum in the denominators are re- 
tained to regulate the endpoint singularity. The masses of the light quarks and the mass 
difference (A) between the h quark and the B meson are neglected. The terms proportional 
to fc^rp, fc|rp in the numerator are dropped, because they are power suppressed compared 
to other terms. In the transverse configuration b-space and by including the Sudakov 
form factors and the threshold resummation effects, we obtain the B ^ P form factors 
as following, 



fi{q^) = ^^T^Cpm^ j dxidx2 J hidhih2dh2ipB{xi-,hi) 

X I [ro (0^(2;2) - 0*(x2)) ■ hi{xi,X2, h, 62) - roXi'qm^cf)'' {X2) ■ h2{xi, X2, h, 62)] 

■asiti)exp [-Sabih)] 
+ [xi {'q(jf{x2) - 2ro0^(x2)) + 4roXi0^(x2)] ■ /ii(x2, Xi, 62, &i) 

■a,(t2)exp[-S„b(t2)]}, (19) 



/2('?^) = 16vrCirm| ^ dx\dx2 J hidhib2dh2ipB{xiihi] 



X 



{X2r] + 1) 0"(x2) + 2ro ( ( ^ - 2:2 ) 0*(a;2) - X20P(x2) + 30" (2:2) 



■hi{xi,X2MM) - roXim% (1 + X2r]) 4>''{x2) ■ /i2(a;i, X2, 61, 62)] ■ as(ti)exp [-Sabiti)] 
+2ro (^ + 1) <P''ix2) ■ hix2,xi,b2,h) ■ as{t2)exp[-Sabit2)]], (20) 



FT{q^) = ^nCprn^B / dxidx2 bidbib2db2{l + rp)'ipB{xi,bi) 



1, (^2 J 
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/ll(Xi,X2,fel,&2) 



■as(ti)exp [-^^{.(ti)] 



+2ror(x2) 



1 + — ) ■ hi{x2,xi,b2,bi] 



as{t2) exp [-S'a6(t2)]|, 



(21) 



where Cp = 4/3 is a color factor, ro = uIq /mB = rn^p/lniBijnq + m^/)], rp = mp/niB 
and rrtp is the mass of the final pseudoscalar meson, rrtq and m'^ is the mass of the 
quarks involved in the final meson. The functions hi and /i2, the scales ti, t2 and the 
Sudakov factors Sab are given in the Appendix A of this paper. One should note that 
fi{q^)-i f2{.(f) and -Ft(q'^) as given in Eqs. f lTUH?!]) do not including the NLO correction. 
In order to include the NLO corrections, the as in Eqs. ( IT9H2T|) should be changed into 
as ■ F{xi, X2, V: f^f: A*) Ci); where the NLO factor F{xi, X2, i]j f^fy l^i Ci) has been defined in 
Eq. dSD. 



B. Semileptonic B and Bs meson decays 

For the charged current B — )■ i^l'vi and i?^ — ?■ K~^l~ui decays, as illustrated in Fig. 1(a) 
and 1(d), the quark level transitions are the b — )■ ul~ui transitions with /~ = (e~, ,t'), 
the effective Hamiltonian for such transitions is 27 



Tieffib ^ uWi) = ^Vub M7/.(l - 75)& ■ h'^il - (22) 

where Gp = 1.16637 x lO^^GeV^^ is the Fermi-coupling constant, V^b is one of the 
Cabbibo-Kobayashi-Maskawa (CKM) matrix elements. The corresponding differential 
decay widths can be written as 



dT{b ^ uWi) _ G^lKfeP - _ /(g2 _ mff j {ml - ml - q^f ^ 



dq^ 1927r3m| (g2)2 y y 4^2 p 

X I (m^ + 2g2) [g2 _ _ ^^)2J ^^2 _ ^ ^^)2J ^2 (^2) 

+3m^ml-mlYF^{q')}, (23) 

where is the mass of the lepton. If the produced lepton is or /i^, the corresponding 
mass terms could be neglected, the above expression then becomes 



^^^^^^^ = ^^A3/2(,2) 1^4,2) |2, 



dg2 192n^m^ 



(24) 



B 



where A(g2) = (m^ + ml — q^Y — 4m^mp is the phase-space factor. 

For those flavor changing neutral current one- loop decay modes, such as i? — )■ PVl^ 
with P = (iiyK) and — )■ K^l~l^ decays, as illustrated in Fig.l, the quark level 
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transitions are the 6 — )■ (s, d)l~l^ transitions, the corresponding effective Hamihonian for 
such transitions is 

10 

(25) 



* i=l 



where q = {d, s), Ci{fi) are the Wilson coefficients and the local operators Oi(/i) are given 

by a 

Oi = {qaCa)v-A{Cfibfj)v-A, O2 = {qaCfi)v~A{cpba)v-A, 

O3 = {qaba)v-A {qWi3)v-A ' = {qJ>p)v-A ^(g^gl)y-A, 



V+A, 



0-j 

09 



^qa^^il + l,)bF,^, 



([7;./) - l6)b] , Oio = ^ {h,7^l) [qY (1 - 75) b] , (26) 

where q = {d, s), q' = {u, d, c, s, b). 

For the decays with b — ?■ sl^l~ transition, for example, the decay amplitude can be 
written as l2j 



^(6^s/+r) = -^^K 



:Vt,{c,o [s7^(l-75)&] [h'l^l] 



+C9^^^(/i) [.-7^(1- 75)6] [[7^/] 



-2m,C7^^^(/i) 



sicr^u^ (1 + 75) & 



(27) 



where Cj^\^) and Cg-^^{^) are the effective Wilson coefficients, defined as 

Cf(^) = C,(/i) + CU^(/x), (28) 
Cfifi) = CM + ^pert(S) + Fld(S). (29) 

Here the term Ypert represents the short distance perturbative contributions and has been 



given in Ref. [29 



Fpert(s) = h{m,, s)Co - h{l, S)(4C3 + + SCs + Cg) 
-^/i(0, s)(C3 + 3C4) + ^(3C3 + + 3C5 + Ce 



(30) 



with Co = Ci + 3C2 + 3C3 + C4 + 3C5 + Ce, s = q^ /m\, rhc = rric/mf, and rhf, = mb/rriB, 
while the functions h{z, s) and /i(0,s) in above equation are of the form 



h{z, s) 



8, mfe 8, 8 4 

In In 2; H 1 — X 

9 fi 9 27 9 



--(2 + x)v/|l-x| 



2 arctan (x = ^ > 11 



/I U, s = m m s H — ztt. 

^ ' ^ 27 9 /i 9 9 



(31) 
(32) 
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The term IldI-s) in Eq. ( I29l) refers to the long-distance contributions from the resonant 



states and will be neglected because they could be excluded by experimental analysis [30 
3l|. The term C'^^^^ in Eq. ( 128|1 is the absorptive part of 6 — )■ S7 and is given by|32[ 



(33) 



with 



Gj{xt) 



Xt (xj — 5xt — 2) 3xj In Xt 



{xt - 1) 



4(xi - 1)^ 



(34) 



where rj = as{mw)/c(s{fJ^) and Xt = ml/m^. 



The differential decay width of 6 — )■ sl'^l is given by[l6|, l33 



dV{h sl+l-) _ GlaLlV^fcPl^/J'v^W) /g' - 1 



512m^7r5 



X 



6m?|CioP(m|-m^)^F^(g^) 



+ (g2 + 2m2)A(g^ 



3g2 



2^7^^^ (m^ - m,)FT(g^ 



+ |Cio|^(g^-4mf) Kq')Fi{q') 



(35) 



where = 1/137 is the fine structure constant. For h — )■ dl^l~ decays, it is easy to 
derive the differential decay width from above equation by a simple replacement Vts — )■ Vtd 
and nis — )■ m^. 

Finally, the effective Hamiltonian for h svv transition is 



b—^suu 



V2 27rsin2(6'vF) 



VtbV,:vxX{xt) [sj^{l - 75)&] [^7m(1 - 75)^] 



(36) 



where 6'wis the Weinberg angle with sin^(6'vK) = 0.231, the function X{xt) can be found 
in Ref. 27(|, while r^x ~ 1 is the QCD factor [27.]. The corresponding differential decay 
width can be written as 



dT{b^suiy) JCl-^'\^X^/^{ml,ml,q^) 



2\l2 



dq' 



(37) 



The factor 3 in above equation arises form the summation over the three neutrino gener- 
ations. For b — )■ duu transition, we can obtain the differential decay widths easily also by 
simple replacements: |V^*| — > l^^^l and — )■ m^. 
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IV. NUMERICAL RESULTS AND DISCUSSIONS 



In the numerical calculations we use the following input parameters (the masses and 
decay constants are all in unit of GeV) js^]: 



A 



(/=4) 
MS 



0.287, = 0.13, /x = 0.16, = 0.21, 

0.236, mB± = 5.2792, rriBO = 5.2795, rriBO = 5.3663, 

1.638 ps, tbo = 1.525 ps, r^o = 1.472 ps, 

0.1396, m^o = 0.135, m^± = 0.4937, m^o = 0.4976, 

1.777, mfe = 4.8, mw = 80.4, = 172. 



For relevant CKM matrix elements we use \Vf 



tb\ 



0.999, \Vf 



ts 



0.0403: 



0.211 ±0.001 ±0.005 |34|. 



hO.OOll 
-0.0007' 



(3^ 



\Vtd/Vt 



ts 



A. Form factors in the pQCD factorization approach 

By using the definitions in Eqs. ( !T6|T7|) and the expressions in Eqs. fll9tl2T]) . we can 
calculate the values of the form factors Fo(g^), F^{q^) and F^lq"^) for given value of 
in the region of < < {Mb — mp)"^. But one should note that the pQCD predictions 
for the considered form factors are reliable only for small or moderate values of q^: say 
< q^ < 12 GeV^. For the form factors in the larger q^ region, one has to make an 
extrapolation for them from the lower region to larger region. 

For the form factor -Fo(q'^) of B/Bg — )■ (vr, K) transition, we make the extrapolation by 
using the pole model parametrization 



Foil') = . (39) 



i^o(O) 

1 — a{q^/rn?B) + h{q^/m\) 

where a, h are the constants to be determined by the fitting procedure. In Table HI we 
list the LO and NLO pQCD predictions for the form factors -Fo(O) and the corresponding 
parametrization constant a and 6 for i? — t- (vr, K) and ^ K transitions extracted 
through the fitting. The first error of (-Fo(0),a, 6) in Table [T] comes from the uncertainty 
of Ub = 0.40 ± 0.04 GeV or ub, = 0.50 ± 0.05 GeV, the second ne is induced by af = 
0.06 ± 0.03 and/or = 0.25 ± 0.15The errors from the uncertainties of 171^' , Vts and 
|Vtrf/Vis| are very small and have been neglected. 

For the form factors and F^lq'^), the pole model parametrization as given in 

Eq. fl39|) does not work, and we have to use other proper parametrization method. In 
this paper, we use the Ball/Zwicky (BZ) parametrization method l^, 35, 36]. It includes 



the essential feature that -F+(g^) and -frlg^) have a pole at = m|*, with B*{1~) is a 
narrow resonance with ms* = 5.325 GeV and ttlb* = 5.415 GeV, which are expected to 
have a distinctive impact on the form factor. 

For the form factors -F+(g^) and FT(g^) of B/Bg — t- (vr, fC) transition, we make the 
extrapolation by using the BZ parametrization 

1 rq^/m^t 

FM') = m) I + ^ I , (40) 

y /^B*^^ 1^1 - q'^/m^J [1 - aq^/m\^ 
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TABLE I. The pQCD predictions for form factors -Fo(O) and the parametrization constant a 
and b for B — )• (vr, K) and Bg ^ K transitions, at the LO and NLO level respectively. The two 
errors come from the uncertainties of uji, or ujb^, and af and/or a^' , respectively. 





Fo{0)lo 




bLO 


B-^-K 


0.22l[!:j]^ ± 0.01 


0.58 ±0.01 ±0.03 


-0.15 ±0.01 ±0.01 


B^ K 


o.27t[!:|];^ ± 0.01 


0.60 ±0.01 ±0.03 


n 1 c;+0.01+0.01 
"•^"-•-0.00-0.02 


Bs^K 


0.22 ±0.03 ±0.01 


o.6i±o.oii[!:[]^ 


-0.16±0.00^[!:[Ji 




Fo{0)nlo 


(INLO 


bNLO 


B^TT 


0.26i[!:[]^ ± 0.02 


0.50 ± 0.01+[];[j^ 


-0.13 ±0.01 ±0.01 


B^ K 


0.31 ±0.04 ±0.02 


0.53 ± 0.0ll[j;[!^ 


-o.i3±o.oit[!:j]? 


Bs^K 


o.26j:[!:[]^ ± 0.02 


0.54 ±0.00 ±0.05 


-0.15 ±0.01 ±0.01 



1.00 




0.00 ^ 







10 



15 



q' (GeV') 



FIG. 2. The pQCD predictions for the form factors Fi{(f') for B ^ it transition, where the 
dots refer to the pQCD predictions for the given points of in the range of < < 12GeV^^ 
while the solid curve stands for the fitted curve at the NLO level. 



where a and r are the shape parameters to be determined by the fitting procedure, the 
same as for the case of Fo(g^). In Table HTl we list the pQCD predictions for the form 
factors -F+(0), -Pr(O) and the corresponding shape parameters (a,r) for B — >■ ij^K) and 
Bg ^ K transitions at the LO and NLO level. In Fig. 2, we show the pQCD predictions 
for the form factors Ft((3'^) for i? — ^ vr transition, where the dots refer to the pQCD 
predictions for each given value of in the large- recoil range of < < 12GeV'^, while 
the solid curve stands for the fitted curve at the NLO level, obtained through fitting by 
using the Eq. (HUI) . 
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TABLE II. The same as in Table [U but for the pQCD predictions for the form factors F+(0), 
-Ft(O) and the corresponding shape parameters a and r at the LO and the NLO level. And the 
two errors comes from the uncertainties of lo^ or cj^sj and and/or a^^ > respectively. 





F+{0)lo 




TLO 


ij — )■ TT 


0.22Iq Q2 ± 0.01 


O.DITq ± 0.01 


0.51 ± 0.00 ± 0.03 


7~> I 

B ^ K 




r\ ().{)() 1 r\ r\'\ 

0.62_Q q]^ ± 0.01 


0.58 ± 0.00 ± 0.03 


Bs ^ K 


[j.ZZ zb U.Uo dz U.Ui 


U.D4 it U.UU it U.Ui 


n p^f^ _i_ n nn 1 +0-04 

U.DD it U.UUlt_g Q3 




Ft{0)lo 




rLO 


B ^ TT 


0.23 ± 0.03 ± 0.01 


r\ r\~\-\) .\]\) i r\ 

0.69lQ;Qj^ ± 0.01 


0.55 ± 0.01 ± 0.03 


B ^ K 


0.30Iq 03 ± 0.01 


O.TIIq q\ ± 0.01 


0.58 ± 0.01 ± 0.03 




(-) nc+0.04 _|_ n ni 
u.^fj Q Q3 HZ u.ux 


n 71 +0-01 + n ni 


Kj.Oiy HI U.UU m u.uo 




-^+(0)ArLO 


OiNLO 


rNLO 


B -^TT 


0.26l[]:[J^ ± 0.02 


0.52 ±0.01 ±0.03 


0.45 ± 0.00^[]:[J^ 


B ^ K 


0.31 ±0.04 ±0.02 


0.54 ± o.oil[J:[!3 


0.50 ± 0.00 ± 0.05 


Bs^K 


0.26l!]:[!3 ± 0.02 


0.57 ±0.01 ±0.02 


0.50 ±0.01 ±0.05 




Ft{0) NLO 


aNLO 


rNLO 


B -^TT 


O.lGtom ± 0.02 


0.65ij]:[!2 ± 0.01 


0.50 ± 0.00 ± 0.00 


B^K 


0.34+°;°^ ± 0.02 


0.67 ±0.01 ±0.01 


0.53 ± 0.00+°;°^ 


Bs^K 


0.28 ±0.04 ±0.02 


0.69 ±0.01 ±0.01 


n c-Q+o.oi+0.04 

'-'•O'J-O.OO-0.02 




FIG. 3. The pQCD predictions for the g^-dependence of Fo.+.tCq^) for i? — )• vr transition, where 
the dots curve and dot-dashed curve shows the LO and NLO part respectively, and the solid 
curve stands for the total value at the NLO level. 



In Figs. ([SllSD, we show the pQCD predictions for the g^-dependence of the form factors 
Fo,+,t(5'^) at the leading order (dots curves) and the next-to-leading order (solid curve) 
for the considered B — )■ (vr, K) and Bg ^ K transitions, respectively. 

From the numerical results as listed in Table [T] and [TTl and the g^-dependence as illus- 
trated in Figs.([3}lSD, 0116 can see that: 
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FIG. 4. The same as in Fig. [3] but for B ^ K transition. 
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FIG. 5. The same as in Fig. [3] but for Bg ^ K transition. 



(i) For the considered B — )■ (tt, K) and Bg ^ K transitions, the NLO pQCD predictions 
for the form factors Fo,+,t(0) agree well with the values estimated from the LCSR 
or other methods 10 . 



(ii) -Fo(O) equals to -F+(0) by definition, but they have different g^-dependence as illus- 
trated by FigsJSHSl We also observe the pattern of the relative strength of the form 
factors: 



<r (0) = ^r^'io) ^ i^cf+ ""(o), 



(41) 
(42) 



which is consistent with the general expectation. 



(iii) The LO part of the form factors dominate the total contribution, the NLO part 
is only around 20%. The form factor -Fo('?^) has a relatively weak g^-dependence, 
but F+(g^) and -Ft((3'^) show a little stronger g^-dependence when compared with 
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B. Decay widths and branching ratios 



By using the relevant formula and the input parameters as defined or given in previous 
sections, it is straightforward to calculate the branching ratios for all the considered 
decays. 

Firstly, in Figs.(|6ll7]), we show the differential decay rates dV /dq^ for the decay modes 
corresponding to the -B — )■ (vr, K) and Bg ^ K transitions. 




q2 (GeV2) q2 (GeV=^) (GeV^) 



FIG. 6. The g^-dependence of the differential decay rates dV/dq^ for the decay processes with 
the S — )• vr and B ^ K transitions. 



By making the numerical integrations over the whole range of g^, we find the numerical 
results for the branching ratios. For the 6 — >■ m charged current processes, the NLO pQCD 
predictions for the decay rates are the following 



|2 



Br{B' -r vr+Z-p,) = (l.42l°:^°(a;,) ± 0.15(a-) ± 0.12^)) x x 10-^ (43) 



IK 



Br{B^ ^ vr+r-p.) = (0.90l°:?^(a;,)l°:°^(4) ± 0.08^)) x x 10-^ (44) 
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10 15 20 

q2 (GeV2) 





q2 (GeV^) 



q2 {GeV2) 



FIG. 7. The g^-dependence of the differential decay rates dV/dq^ for the decay processes with 
the ^ K transitions. 



Br{B- n^rui) 



Br{B- TT^T-Ur] 



Br{Bl K+T-Vr) 



(7.63l?i?l 



+0.82/ 7rN+0.66, 
-0.78 1"2 J -0.63' 



o)) 



IK 



ub 



(l.27+o;26(wfes 



+0.47, 
-0.41* 



|0.0038|2 

\Vub? 

I0.0038P 



X 10" 



X 10" 



+0.14/ X\+0.10, 



-0.13 



a. 



-0.09 



where the first error comes from the uncertainties of Wf, 



+Q.m( „K\+0.&2, 



-0.66 



a. 



-0.59 



IK 



10.00381 



IK 



10.00381 



X 10" 



X 10" 



(45) 
(46) 
(47) 
(48) 



the second error are induced by the variations of 
the third error comes from the uncertainties of itiq 



K 



0.06 ±0.03, 
1.4±0.1G'eK 



0.40±0.04 or = 0.50±0.05, 
0.25±0.15and 
:1.6±0.1GeK. 

For B^ — )■ 7T^l~ui and B~ — )■ ir^l'vi decay mode, their branching ratios have been 
well measured by BaBar, Belle and CLEO Collaborations [l|-|3|. The new BaBar 
measurement [l]] and the new world average [s^l are the following 



Br (5° ^ Ti^rui) 



(1.41 ± 0.05(s?/st) ± 0.07(stat)) x BaBar[l], , . 

(1.44 ± 0.05) X 10-^ PDG2012[37], ^ ' 

Br{B- n^ri^i) = (7.78 ± 0.28) x 10"^ PDG2012[37], (50) 

On the other hand, we know that one can extract out the magnitude of the CKM matrix 
element Kfe by comparing the theoretical prediction for Br{B — )■ Trlu) with the data. 
Based on the measured partial branching fraction for B — )■ vr/z/ in the range of < 
< 12GeV^ and the most recent QCD light-cone sum-rule calculation of the form factor 
F+(g^)j38[, BaBar Collaboration found the result |li 



iKfel = (3.78±0.13(exp. 



)t°oil{theor.)) X 10" 



(51) 
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where the two errors refer to the experimental and theoretical uncertainties. 

From the differential decay rate as given in Eq.f l2^ and the pQCD calculation of the 
form factor Fj^{q^) at the NLO level, we make the numerical integration over the whole 
range of < < {Mb — Tn^^Y, compare the obtained branching ratio for — > n^l^ui 
with the new BaBar measurement as given in Eq. (H9|) . and derive out our estimation for 
the magnitude of Vub- 

\Vub\ = (3.80l|]i°(u;,) ±0.20(aJ)l°:{^K)) x lO'^ = (3.80l°:^^(tW.)) x 10-^ (52) 

It is easy to see that our estimation for the central value and the uncertainty of \Vub\ 
agrees very well with the BaBar result as given in Eq. ( ISTl) . 

For other neutral current processes, after making the numerical integration over the 
whole range of < < (M^— m^)^, we find the NLO pQCD predictions for the branching 
ratios. The pQCD predictions at the NLO level and currently available data are all listed 
in Table lllli The first error of the pQCD predictions comes from the uncertainties of 
Ub or ujbs, the second one from af^ and/or a^'^ , and the third one is induced by the 
uncertainties of the chiral mass ml'^ . 



TABLE in. The NLO pQCD predictions for the branching ratios of the considered decays 
with / = (e, (i) and currently available experimental measurements |39| and the world 

averages " 



371 ]. The upper limits are all given at the 90%C.L. 



Decay modes 



Br{B^ 



Br{B- 
Br{B' 
Br{B- 



Br{B^ 
Br{B^ 
Br{B^ 



Br{B- 
Br{B- 
Br{B- 



Br{B^, 
Br{B^^ 
Br{B^ 



vrOr+r- 



RH+l-) 



K-in-\ 

K-T+T- 

K-vv) 



K'^l+l-) 



NLO pQCD predictions 



o.9i: 

0.28: 
7.30 



±0.10 ±0.05 



-0.26 
-0.19 

± 0.02 ± 0.03) 



X 10" 
X 10" 



-2.07H 
-1.54- 



0.79H 
0.74- 



0.63^ 
0.61J 



X 10" 



1.95: 
o.6o: 

1.57 



-0.55+0.21+0.17 
-0.41-0.20-0.16 
-0.16+0.04 
-0.12-0.03 
-0.44+0.17+0.14 
-0.33-0.16-0.13 



X 10 

± 0.06) X 10 

10-7 



-8 



5.1 



1.20" 



+1.5+0.5+0.4\ 



X 10 



-7 



-1.1-0.5-0.4 
0.32+0. 07+0. ll^ 
0.25-0.07-O.loy 
.2+0.4+0.3^ ^ -in-6 
.9-0.4-0.3' ^ 



X 10" 



5.50: 
1.29: 
4.42 



-1.59+0. 57+0. 42\ 
-l.18-0.55-0.4i; 
-0.35 
-0.26 

-1.28+0.46+0.34\ 
-0.95-0.44-0.33^ ^ 



X lO-'^ 
±0.08 ±0.11) X 10- 
10-6 



1.633 

0.43 

1.31 



+0.54+0.18+0.12- 
-0.38-0.17-0.12, 
0.13+0.03+0.04^1 
0.10-0.03-0.04^ 
0.43+0.14+0.10^ 
0.31-0.13-O.lOy' 



I X 10" 

X 10-^ 
X 10-' 



Data 



< 1.2 X 10" 



< 2.2 X 10"^ 



< 4.9 X 10" 



< 1.0 X 10-^ 



(4.7^°:^) X 10 



< 5.6 X 10" 



-7 



(5.1 ±0.5) X lO-'^ 



< 1.3 X 10" 



From the NLO pQCD predictions for the branching ratios of all considered semi- 
leptonic decays of B and B^ meson, as listed in Eqs. (I44ti47p and Table llllt we have the 
following points: 

(i) The branching ratios of the charged current processes B — > ttIv and Bs — >■ Klv 
are all at the order of 10"^. For B^ — )■ t[^1~vi and B~ — )■ tt^I^ui decay modes, the 
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pQCD predictions for its branching ratios as shown in Eqs. ( 143)1451) agree very well 
with the data as given in Eqs. f HUllSU]) . For other charged current decay modes, the 
pQCD predictions as given in Eqs. fl44f461148|) will be tested by the LHCb and the 
forthcoming Super-B experiments. 

For the neutral current 5° K°l+l- and R-l+l' decays, the NLO pQCD 

predictions for their branching ratios agree very well with currently available exper- 
imental measurements. For other neutral current decays, the NLO pQCD predic- 
tions are all consistent with currently available experimental upper limits and will 
be tested by LHCb and the fothcoming Super-B experiments. 



(iii) Because of the strong suppression of the CKM factor iVtrf/Vt^p = 0.211^ |3J], the 
branching ratios for the decays with 6 — )■ d transitions are much smaller than those 
decays with the 6 — )■ s transitions. Furthermore, the branching ratios of — ?■ Pvv 
are almost an order larger than their corresponding decay modes — t- Pl^l^ 
partially due to the generation factor Ng = 3. In order to reduce the theoretical 
uncertainty of the pQCD predictions, we defined several ratios R^, Rc and Rni,N2,N3 
among the branching ratios of the considered decay modes. 

The NLO pQCD prediction for the ratio Ri, is of the form 



Ru 



Br{B^ ttH+I-) Br{B^ KH+1-) 



Br{B- n-l+l-) Br{B- R-l+t 
BriB^^ K^uu) 



for / = {e,fi). These relations will be tested by experiments. 



(53) 



(iv) Because of the large mass of r lepton, we found that the considered B/Bg decays 
involving one or two r's in the final state have a smaller decay rates than those with- 
out r. The pQCD predictions for the ratios Rc and Rni,N2,N3 of the corresponding 
branching ratios of relevant decays are the following 

5r(i?°) ^ P+l-i^i) , , 

Br{B^ n^l+l-) Br{B- n-l+l-) 
" BriB^ ^ vrV+r-) Br{B- ^ n-r+r-) ^^^^ 

Br{B^ KH+1-) Br{B- R-l+l-) 
~ Br{B^ ^ K^T+r-) Br{B- ^ R-t+t-) ^^^^ 
_ Br{Bl-^KH+£l ^ 

- Br{BO ^ ROr+r-) ^ ^^^^ 

for I = (e, fi). These relations will be tested by LHCb and the forthcoming Super-B 
experiments. 
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V. SUMMARY AND CONCLUSIONS 



In this paper we calculated the branching ratios of the semileptonic decays B — )■ 
{n, K){l~^l~ JUjUu) and Bs — t- K{l^l^ ,lu,uu) in the pQCD factorization approach. We 
firstly evaluate the B — )■ {7t,K) and Bg K transition form factors -Fo^+,t(q'^) by em- 
ploying the pQCD factorization approach with the inclusion of the next-to-leading-order 
corrections, and then we calculate the branching ratios for all considered semileptonic 
decays. Based on the numerical results and the phenomenological analysis, we found the 
following points: 

(i) For the B — )■ (tTjK) and Bs ^ K transition form factors Fo^+,t(q'^), the NLO 
pQCD predictions for the values of Fo,+,t(0) and their g^-dependence agree well 
with those obtained from the LCSR or other methods. The NLO part of the form 
factors in the pQCD factorization approach is only around 20% of the total value. 

(ii) For the charged current B^ — t- n^l^ui and B^ — t- n^l^ui decays and the neutral 
current B^ — )■ K^l^l~ and B~ — )■ K^l^l^ decays, the NLO pQCD predictions for 
their branching ratios agree very well with the measured values. 

(iii) By comparing the pQCD predictions for Br{B^ tt^I^iji) with the mea- 
sured decay rate we extract out the magnitude of the CKM element Vub- \Vub\ = 
{3.80+'^ofo{theor.)) x lO'^. 

(iv) We also defined several ratios of the branching ratios Ru, Rc and Rni,N2,N3, 
and presented the corresponding pQCD predictions, which will be tested by LHCb 
and the forthcoming Super-B experiments. 
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Appendix A: Related functions defined in the text 



In this Appendix, we present the functions needed in the pQCD calculation. The 
threshold resummation factors St{x) is adopt from Ref. 13l |: 



2i+^T(3/2 + c) , 



(Al) 



and we here set the parameter c = 0.4. The hard functions hi and /i2 come form the 
Fourier transform and can be written as 40 



+9{b2 - hijlo (i/x^ms6i) Kq {y/x^mBh2) St{x2), (A2) 
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h2{Xi,X2MM) 



biKi {^yxlX21]mBbl) 



9{bi -b2)Io (a/^"^b&2) Kq {y/x^niBbi 



2y/x^x^mB 

+9{b2 - bi)Io {^/x^mBbl) Kq (v'a^ms&2) St{x2), 



(A3) 



where Jq is the Bessel function and Kq, Ki, Iq are modified Bessel functions. 

The factor exp[—Sab{t)] contains the Sudakov logarithmic corrections and the renor- 
mahzation group evolution effects of both the wave functions and the hard scattering 
amplitude with Sab(t) = SB{t) + Sp(t), where 



Ssit) = s ( 

Sp{t) = S ( X2^,b2 



+ 2 



l/b2 



(A4) 
(A5) 



with the quark anomalous dimension 7^ = —as/ir. The functions s{Q,b) are defined by 
l3 



s{Q,b) 



Ad) 



2(3, 



-gin ( 1 



2(3 



,4(2) ^(1) 



In 



4/3? Kb 
In (I 



' 4/3? 



ln(2g) + 1 ln(26) + 1 



+ 



In2(2g) -ln2(26) , (A6) 



where the variables are defined by g = ln[Q/(V2A)], b = ln[l/(6A)], and the coefficients 

A^'^ and (3i are 



/3i 



33 -2n 



/ 



12 



/32 



153 - 19n 



/ 



24 



3' 



^(2) = ^ + ln(e^^/2) 



9 



27 



(A7) 



here, is the number of the quark fiavors, the 'Je is the Euler constant. The hard scales 
ti in the equations of this work are chosen as the largest scale of the virtuality of the 
internal particles in the hard 6-quark decay diagram, 



ti = max{y/x^mB, I/&1, 1/&2}, h = msix{y/xrr]'mB, I/&1, 1/&2} 



(A8) 
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